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Abstract 

A continuous version of the hierarchical spherical model at dimension d = 4 is investigated. 
Two limit distribution of the block spin variable X 7 , normalized with exponents 7 = d + 2 
and 7 = d at and above the critical temperature, are established. These results are proven by 
solving certain evolution equations corresponding to the renormalization group (RG) trans- 
formation of the O(N) hierarchical spin model of block size L d in the limit L J, 1 and N — > 00. 
Starting far away from the stationary Gaussian fixed point the trajectories of these dynam- 
ical system pass through two different regimes with distinguishable crossover behavior. An 
interpretation of this trajectories is given by the geometric theory of functions which describe 
precisely the motion of the Lee— Yang zeroes. The large-iV limit of RG transformation with 
L d fixed equal to 2, at the criticality, has recently been investigated in both weak and strong 
(coupling) regimes by Watanabe [Wj. Although our analysis deals only with N = 00 case, it 
complements various aspects of that work. 



1 Introduction and Statement of Results 

We continue the investigation starting in |CMj . In the present work we give a geometric interpreta- 
tion to certain trajectories of a first order partial differential equation related to the renormalization 
group transformation (RGT) of a d-dimensional hierarchical spherical model. 

*Partially supported by CNPq and FAPESP. E-mail: marchett@if.usp.br 

T Supported by FAPESP under grant #98/10745- 1. E-mail: william@if.usp.br. 
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Motivation. The hierarchical O(N) spin model, with L d = 2 sites per block, has been recently 
studied by renormalization group in both weak and strong regimes by Watanabe [W]. Starting 
from the uniform "a priori" measure supported in the iV-dimensional sphere of radius y/N, the 
critical trajectory of the RGT has shown to converge to the Gaussian fixed point for sufficiently 
large N. To control such trajectory, which starts far away from the fixed point, the exactly solved 
0(oo) trajectory has been used together with two key ingredients: reflection positivity and the Lee- 
Yang property of single-site "a priori" measures. The former ingredient gives uniform convergence 
of 0(N) trajectories to O(oo) trajectories. The latter property has been previously employed by 
Kozitsky [K] to establish two central limit theorems. Watanabe's analysis, based in his joint work 
with Hara and Hattori [HHW] on the critical trajectory for the hierarchical Ising model (N = 1), 
in contradistinction to Kozitsky's, and most of the previous studies of this model, does not restrict 
the space of "a priori" measures to a neighborhood of the Gaussian fixed point and is able to deal 
with the borderline d — 4 case. 

Although the analysis of the RGT with L d > 2 fixed is expected to be simplified considerably 
in the L J, 1 limit (see e.g. [F]). none of the above mentioned results can be carried to the 
limit as the two key ingredients do not hold if L d is not an integer. In order to establish, in the 
local potential approximation (L [ 1), a weak convergence of the hierarchical O(N) Heisenberg 
equilibrium measure to the corresponding spherical equilibrium measure as N —>■ oo an entirely 
new method of analysis has to be developed from scratch. 

In the present investigation we establish central limit theorems for the four-dimensional hierar- 
chical spherical (N = oo) model at and above the critical temperature. Our results are achieved in 
the local potential approximation that reduces the renormalization group equation to a nonlinear 
first order partial differential equation. A geometric function interpretation of the O(oo) trajectory 
is thus given with the help of an explicit solution obtained by the method of characteristics. It 
follows from our analysis that the Lee- Yang zeroes reach a limit distribution as the Gaussian fixed 
point approaches but their support moves away to infinity. 

The model. The hierarchical Heisenberg model on a finite box A K = {0, 1, . . . , L K - l} C Z d of 
size n = L dK is given by the O(N) invariant equilibrium measure 

dvl N \x) = -^exp {I (x,Ax)Af[dai N) ( Xj ) (1.1) 

V ) j = i 

where ) denotes an element of the configuration space Q n = M. N x • ■ ■ x R N ; A = J®I 

the tensor product of the hierarchical coupling matrix J (whose quadratic form 

(s,~Js) Ak = -(L-l) f>- 2fc £ {B k s)l , 

k=l reA K ~k 

ie{o,...,i-i} d 

coincides with Dyson's hierarchical energy [D] when there are L d = 2 sites per block)H with the 

1 The factor L — 1 is chosen so that the hierarchical Laplacean converges, as L J, 1, to a continuum hierarchical 
Laplacean (see [F1 [CM] V 
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N x N identity matrix /; er (x) the "a priori" uniform measure on the iV-dimensional sphere 
\x\ 2 = j3N of radius y/JJN with (5 the inverse temperature. 

Recursion relations. The invariance of J under block transformation (jl.2p allows to establish a 
recursion relation: 

a W( x ) = J_ e ^L-l)^/2(N) + _ # a (N) ( r7/2 



i(L^x) (1.3) 



L d — times 



on the space of single-site "a priori" measures in R with initial data <Tq (x). Here, * denotes the 
convolution product 

p * T](x) = / p(x — x') drj(x') , 



Ck is chosen so that crj^ is a probability measure and 

/ 1 if 7 = + 2 

The "a priori" measure at the step k is defined by integrating (11 .11) over f2 n with the value of 
fc-th block spin fixed: 



1 5{{B k ®l)y-x) dui N \y) = -^ex P \Ux,A X ) UL \ ]J da[ 

J "riW-H V ) „_1 



L d(Jf-fc) 

(AT) 



J L d(K-k) «• ' j = l 

is a marginal measure on fl^K--*) that preserves the form (II. ip . 
In terms of their characteristic functions 

equation (ll.3p reads 

1_ L + l A \ f ,(m A Ld 

4 

for /c > 1 with 



^exp (1.6) 



*o™(*) = , ^^SZT ^-r (v^kl) - rf°(W) ■ (1.7) 
(#V|z|/2) W/i V ) 

Here, exp (£A) is the semi-group generated by the iV-dimensional Laplacean operator A = d 2 /dzf + 
• • • + d 2 /dz%, Nk is chosen so that 0^(0) = 1 holds for all k = 1, . . . , K and J a {x) is the Bessel 
function of order a (see eq. (20) in Chapter VII of |CHj for an appropriate integral representation). 
Note that (jy^\z) = ^4 (r) depends only on r = \z\ = yjz • z. 

Thermodynamical functions. The macroscopic behavior of the model is described by the limit 
distribution of the block variable 

1 n 



where 7 is chosen in order the limit law to be attained. The characteristic function associated with 
the block variable X^ N with 7 = d + 2 is given by 



x 



= J exp (ix ■ z) da { p{x) = (p]p(\z\) . 

As n goes to infinite, N converges in distribution to Xjj if <f jP(r) converges at every point 
r > to a function (p( N \r) that is continuous at r = 0, by continuity theorem (see e.g. [D]). The 
convergence of = lim^oo to the equilibrium measure v of the spherical model is more 
subtle and we analogously employ: Xjj is said to converges to X 7 in distribution if 

lim (ipW(y/Nr)) 1/N = ip(r) 

exist for every point r > 0, is continuous at r = and coincides with the corresponding characteristic 
function of the spherical model. The reescaling is seen to be necessary already at the initial function 
fo ( r ) ( see Proposition 12. II) . 

The statements about convergence are independently of which order both limits n — > 00 and 
N — > 00 are taken. This has been shown in [CM] adapting a method employed by Kac and 
Thompson [KT] for the hierarchical equilibrium measure fim with 7 = d, L d > 2 an integer and 
(3 different from the critical inverse temperature f3 c = /3 c (d, L) of the hierarchical spherical model. 
In |CM] . X d is shown (see Theorem 2.3 and Remarks 4.2) to be Gaussian with mean zero and 
variance where \i = fi(j3) is implicitly defined by 

0= f -r-- — dg(X) (1.9) 
J 

where q{X) is the density of eigenvalues (counted multiplicities) of the hierarchical Laplacean 

L - 1 , 

—Ah = 1 — J ■ Note J is not invariant under translation by a vector in Z , property 

L 2 — 1 

that is required for coupling matrices in [KTj . Some statements about hierarchical spherical model 
hold also in the limit as L I 1, in which case (I1.9P reads (see Section 3 of [CMj lR 

1 - T = -2A*ln f 1 - J (1.10) 



4 r \ 2/2 / 

for d = 4. Central limit theorems are established in the present work directly from the L J, 1 limit. 
Local Potential Approximation. Let 

[/(*,*) = - In tf^) (1.11) 

be defined for £ = In L. As — >• 00 together with L J, 1 so that k In L is kept fixed at a positive 
real number t, (11.41) converges to 



c 7 (t) 



1 if 7 = d + 2 
e -2 * if 7 = d 



2 By monotonicity, there exist a unique solution = < denned for < (3 < 4. 
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and we have 

U(t, z)-U(t-]nL,z) 



U t = lim 

Z41 InL 



lim — •{ - In 

k— >oo t 



p td/k' 



Jexpf-^Aj^He-^))"' 
Consequently, (11. lip satisfies the initial value problem 

tf t = (At/ - \U Z \ 2 ) +dU-lz.U, + ^AU(t, 0) (1.12) 



with 



U(0,z) = -\n^>(z) . (1.13) 

The last term in the right hand side ensures that U(t, 0) = for all t > 0. Note that this property 
is satisfied by the initial condition because of the normalization J a^ N \dx) = (J)q N \o) = 1. 

We shall prove two limit theorems (Theorems 12.21 and 14.11) summarized as 



1 f I I 2 

lim lim —U(t,VNz) = \ |2 



|z| 2 /2/i if (3<f3 c 



uniformly in compact subsets of ( G C with 9?e (£) = — |^| 2 . The first, when the sum (II. 8ft is 
normalized with abnormal exponent j/d — 1 + 2/d, holds at the critical point 

2d 

P = &W = J=2> (1-14) 

d > 4. The second, for normal exponent 7/d = 1, holds for any /3 < /3 c (c?) and d > 2. In both cases 
only the borderline d = 4 will be considered for brevity. 

Conformal mapping. Although continuity at |z| = suffices for these limit theorems, the 
"characteristic function" lini7v-*oo ex P (^j^U^jy/Nz)^ is shown to be an analytic function that 

converges, as t —* 00, to an entire function. In addition, thanks to an explicit solution of the initial 
value problem (11.121) and (11.131) at iV = 00, the whole trajectory can be described by the geometric 
function theory. 

The initial value (11.131) is a function of \z\ 2 and equation (11.12p preserves this property. So, we 
define 

u{t,x)= lim —U{t,y/Nz) (1.15) 

for x = — \z\ 2 and let, for each t > 0, the partial derivative u x (t, Q of u be extended as an analytic 
function of £ = x + iy with y > 0. We prove in Theorem 13.21 that u x (t, £), t > 0, map the upper 
half-plane H conformally into a decreasing family of open convex sets 

u x (t,M) = Q t cQ = u x (0,H) 

contained in H, and there is a one-to-one and onto relation between this family and the trajectory 
at the critical inverse temperature (3 C {4) = 4 converging to the Gaussian fixed point. Analogous 
theorem holds for the trajectory corresponding to normal fluctuations. 
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The boundary of Q t is the union of a segment I a := [— a, 0] extending from a point —a = 
—a(t) < up to the origin over the real line and a convex curve q = h(t,p), p G I a , with 
h(t,—a) = h(t,0) = 0. h(t,I a ) = {h(t,p),p G I a } encodes all informations about 6 since it 
corresponds to the image of a branching cut of u x (t, Q. The principal branch of u x (t, () belongs to 
the Pick class of analytical function and admits to be represented as 

M*,0 = -i + £ - «W*,A) (Lie) 

where dp = p dX is absolutely continuous (with respect to Lebesgue) Borel measure. Although 
f ll.l6p is not a canonical representation, 

p (t, A) = — lim S (u x (t, A + it])) 

7T 77J.O 

holds as well. Denoting by £(t) = (— 00, —d(t)) the support of p in (jl.lfip . we have 

-a(t) = u x (t, -d(t)) 



h(t,I a ) = S (u x (t, E(t) + i0)) 



The support of p{t,X) determines the location of the Lee- Yang zeroes as it can be seen 
by representing tp^\r) into a infinite canonical product (for <fo*\r), see proof of Proposition 12. II) . 
By (11 .lip and (11. 15ft . these zeroes are poles of u x (t,Q that become dense over the semi-line £(t) 
as — > 00. As t goes to 00, a(t) — > 3/2, d(t) — > cxd leading S(t) to an empty set as all Lee- Yang 
singularities are expelled to infinite. As a consequence, u x (t,C) — ^ — 1 uniformly in each compact 
set of C. 

The motion of the Lee- Yang zeroes can be attained from the moments of their distribution 
[Nj . The moments satisfy an infinite system of ordinary first-order differential equations which is 
reduced in [HHWt IW] to a finite system by Lee- Yang inequalities. The presence of one-dimension 
unstable manifold makes this system very sensitive to truncation and no simplification occurs in 
the limit iV — > 00. This has to be contrasted with the simple geometric analysis in Section [3] from 
which the dynamics of Lee- Yang zeroes can be described globally. 

Outline. In Sections [2] and H] we prove Theorems 12.21 and 14.11 which are Gaussian limit laws for 
the spherical model on the local potential approximation. Section [3] presents an interpretation of 
explicit solution of the associate nonlinear first order partial differential equation according to the 
geometric function theory. A conclusion with final remarks is given in Section [5j 



2 Central Limit Theorem 



The radial equation. The initial value (I1.13P is a function of \z\ 2 = r 2 and the spherical sym- 
metry is preserved by the evolution equation (I1.12p . So, it suffices to take into account the radial 
component of z ■ d/dz and A, respectively given by rd/dr and 

»0"»A 9 
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Defining 

u^(t,x) = ^U(t,^Nz) (2.1) 



for x = — | z | , the initial value problem (11.121) and (11.131) for 7 = d + 2 reads 

u[ N) = ^*u2? + «f > - 2x (n^) 2 - lxu{ N ) + eft*™ - «£">(*, 0) (2.2) 
with (0,x) = t/(0, y/Nz)/N. As iV -> 00 , the initial function converges to a limit: 
Proposition 2.1 

lim «W (0, X ) = f Z^-r^ dx' ee u (x) (2.3) 
and t/ie convergence is uniform in any compact set of the slit plane C\(— 00, —1/4/3]. 

Proposition 12. II is proven in Section [3) Watanabe established (12.31) writing Uq(x) as a continued 
fraction of Gauss (see Lemma 4.1. of [ Waj ) . Additional properties are obtained by taking into 
account that u' is an analytic function of the Pick class Pi{p) which is able to be continued across 
the interval I((3) = (—1/4/3, 00). 

Viscosity limit equation^. Taking iV — >• 00 in f |2.2j) we are led to a first order partial differential 
equation for (11.151) 

Ut = u x — 2xu 2 x — ^xUx + du — u x (t, 0) (2.4) 

which can be solved by the method of characteristics. To avoid dealing with a nonlinear equation 
we apply the Legendre transformation to (12.41) . Let 

w(t,p) = max (xp — u(t, x)) = xp — u(t, x) (2.5) 

x>0 

be the Legendre transform of u with respect to x where x = x(t,p) is attained at the value x for 
which 

p = u x (t,x) (2.6) 

has a solution for every t > and p in a certain domain depending on t. 

Assuming w(t,p) continuously differentiate and uniformly convex function of p such that 
linip^oo w(t,p)/ \p\ = 00 holds for all t > 0, the original function u(t,x) can be recovered by 
inverse Legendre transformation 

u(t, x) = max (xp — w(t,p)) = xp — w(t,p) (2.7) 



where p = p(t, x) solves x = w p (t,p) for p. Note that, by differentiating (12.51) with respect to t and 
p together with (12. 6p . we have 



w t = -u t 



w 



V 



x + (p — u x (t,x)) Xp = x . (2.8) 



3 1/N plays the role of viscosity since it is in front of the Laplacean as in the hydrodynamic equation of incom- 
pressible fluid. Viscosity solution (or limit) also refers to a method for obtaining "weak solutions" of semilinear first 
order partial differential equations (see e.g. [E]). 
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Hence, w p solves equation (12.61) for x. We are going to show that w p (t,p) is a monotone increasing 
function of p for every t > therefore, w(t,p) is convex and a well defined Legendre transform of u 
which, by (12.71) . is also uniformly convex. It follows by duality of the Legendre transformation that 

p(t,x)=u x (t,x) (2.9) 

which, in view of the presence of u x (t, 0) in (12 .411 . gives 

PX 

u(t,x) = / p{t,x') dx' . (2.10) 



JO 

Using 7 = d + 2 together with (I2.5P and (12.81) . equation (12.41) becomes 

w t = -p + 2p (1 + p) w p + dw + po 

where p = po(t) is implicitly defined by the equation = w p (t,p). Writing v = w p = x we arrive, 
by differentiating both sides of the above equation with respect to p, at the following initial value 
problem 

v t -2p(l+p)v p = -l + (j + ip)v (2.11) 

with 

v{0,p) = ± + -£ i = v o (p). (2.12) 
Note that v (0,p) = x(0,p) is the value x that solves (12. 61) at t — 0: 

p = u x (0, x) = u' (x) = ~^=r ( 2 - 13 ) 

1 + vl + Apx 

by ([Z2D and (I23|) . 

Main result. Our main result of this section is as follows 

Theorem 2.2 Equations h2.11\) and k2. lty) with d = 4 (j — 6) are solved by 

v(t, P ) = ^ + - 2 -^-e*- l -±^ln(l+p-pe 2t ) . (2.14) 
2p p z Ap z p A 

At (3 = (3 C = 4, there is a unique solution p = p(t, x) of 

v(t,p)=x (2.15) 

holomorphic in a neighborhood of origin, that converges, as t —* oo, to —1 in every compact set of 
C. Together with equations Ii2.10\) and \2. this implies convergence to the Gaussian equilibrium 
solution of M.12\) : 

lim lim — U(t,\/~Nz) = \z\ 2 

t->ooN-*oo N 

uniformly in compacts. 
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Remark 2.3 The use of Legendre transform in the renormalization group transformation for the 
0(N) Heisenberg model in the large-N limit goes back to Shang-Keng Ma's work (see /Ma]/ and 
references therein). It is also reminiscent of the method of Laplace (see eqs. (3. 1.13)- (3. 1.16) 
of IKKPStf ). In ref. IWf . Watanabe solved the discrete flow equation U.6\) with L d = 2 in the 
N — > oo limit and partial differential equation is employed only for the heat semigroup part in 
U.6\) . Theorem \2.2\ extends Watanabe's result to the flow equation hi. 12) at the L [ 1 limit. 



Remark 2.4 Theorem \2.2\ treats the border case d = 4 but holds for any d > 4. The proof of 
the theorem can also be adapted to deal with the convergence to nontrivial equilibrium solutions of 
firm) at 13 = (3 c {d), given by (tt$ , for 2 < d < 4. 



Proof. Theorem 12.21 will be proven by solving (12.1 ip along the characteristics p(t) = p(t; po) (see 
e.g. [E]). Writting V(t) = v(t,p(t)), equation ( 12.111) is reduced to a pair of ordinary differential 
equations 

p = —2p(l+p) 

V = -l + (6 + 4p)V (2.16) 
satisfying initial conditions p(0) = po and 

V(0) = Vo = v (po) . (2.17) 

Integrating the first equation of (12.161) 

„ P'(l + P') Jpo VP' 1+P7 Jo 
gives 

— 2t 

P(t) = ~ t ■ (2.18) 

The second equation of (12.161) is a nonhomogeneous linear equation. The homogeneous equation 
V = (6 + Ap) V can be integrated: 

V(t) = V exp ^6t + 4 J p{s) ds 

= V e 6t (1 + Po - Po e^) 2 . 

Using the variation of constants formula (see Theorem 3.1 of jCL] ). the solution to the second 
equation of (12 . 161) is given by 

V(t) = e 6t (1 +p - Po e- 2t ) 2 (Vo - J ) (2.19) 

with 

Jo = 



e" 6s ds 



,1 +Po ~Poe~ 
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by changing variable £ = e , given by 

1 



;i + Po) 2 -— r + 2(1 + po) In (1 + po - PoC) + PoC 

1 + Po - PoC 



exp(-2t) 



After some manipulations together with (12.171) and ( 12. 12jl . this gives 



V - J = ^ - -L + () + Q + l + g£ln(l + Po -p„e-") . (2.20) 

4p 2 2pg 2pg(l+p -p e- 2 *) 2p 2 pg v 



Equation (12.141) follows by plugging this result into (12. 19[) with p as a function of i and p: 



2/ 

pe 2 ' 

Po(^p) 



1 + p — pe 2t 
obtained by solving (12.181) for p . 

The inverse function theorem. We now solve equation (12.151) for p at the critical point f3 = 
P c (4) = 4. By (I2.14p . it can be written as 

xp 2 - - - 1 = _i±^ l n (i +p - pe 2t) = g (t p ) . (2.21) 
2 p ' 

The first of two ingredients we need is 

Lemma 2.5 For p < (e 2t — 1) 1 , g is a monotone increasing function of p diverging to — oo 
logarithmically as p — > — oo and satisfying g(t, —1) = and g(t, 0) = (e 2 * — 1). 

Proof of lemma. Clearly, g is well defined function of p for 1 + p — pe 2t = 1 — p (e 2a * — 1) > 
with logarithmic divergence at p = — oo. We have, by an explicit computation, 

where 

/(to) = ln (1 - w) + — l = h(w)-l. (2.22) 

1 — w 

If f(w) > for all w < 1 then g p (t,p) > in the domain p < (e 2 * — 1) 1 and the monotonicity 
statement is proven. In fact, h(0) = 1 and 

h\w) 



l + wf 



implies that w = is the absolute minimum of h proving an equivalent statement: h(w) > 1 for 
w < 1 different from 0. 

□ 

For x < 0, the quadratic polynomial 

Q(x,p) := xp 2 - - - 1 
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in the left hand side of (12.211) is bounded from above by a linear function: 

Q(x, P )<Q(0,p) = -^-l, 

and attains its maximum value 1 at p max = — . Since p max — > as x — > — oo, there is a value 

16x A x 

^max = ^max(^) such that no real solutions of (12.1 5p exist for x < x max . On the other hand, as the 
graph of g(t,p) intercepts the graph of Q(x,p) in two points (one point) for any > x > x max (t) 
(x > 0) and t > 0, there exist at least one real solution of (12.151) for x > s max (see Figure ). We 
shall discard the solution associated with the second point of interception since it diverges at x — 0. 

Now, let t > and let x and p be real parts of numbers in C: z = x + iy and 77 = p + iq. 
Although the solution 77 = r](t, z) of z = v(t, 77) is a multivalued function of z, only one branch, 
denoted by fj(t,z), is regular at z — 0. Note that fj(t, 0) exists for all t > and is a real valued 
monotone increasing function of t > satisfying —2 < f](t,0) < — 1 as the graph of g(t,p) always 
intercepts the straight line Q(0,p) = —p/2 — 1 at some negative point p*(t) within that range (see 
Figured]) and p*(t) = fj(t, 0) by definition. 
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Figure 1: Graphic of intersections of g(t,p) (solid line) with Q(x,p) for x = (dotted line) and 
1 / (dashed lines) 

It follows that v (t, rf) is holomorphic in 3?e (77) < with v (t, fj(t, 0)) = and 

11 l/l\e 2 '-l 
Mt,*t,0)) = ^ + - (e 2 ' - l)fl) +^ (l + - J ; _ (e2 , _ t) - > (2.23) 

and these are the assumptions of our second ingredient (see Theorem 9.4.1 of [Hi] for a proof) 

Theorem 2.6 Let R > r > and fj G C be such that v(t,rj) is holomorphic in Dr(j}) = 
{rj G C : 1 77 — 77 1 < R}, v(t, fj) = 0, v v (t, fj) > and v(t, rj) 7^ for < 1 77 — 77 1 < r. T/ien i/ie 
contour integral 

:= ^i v ^M^-z dr] 
11 



where c io = {i] G C : \q — fj\ = p} for some p < r, defines a holomorphic function in {z : \z\ < m} 
where 

m = mm \v(t, fj + pe l6 ) | . 
Moreover, 7] = fj(t,z) is the unique solution of z = v(t,r]) regular at z = in this domain. 

For fixed t, let R = R(t) be such that Dn(fj(t, 0)) C {3?e (77) < 0} and note that we can 
always take R large enough to include 77 = —1. Let r < R be so that v(t,rj) 7^ for < 
1 77 — 77(2, 0)| < r. This is always possible by continuity in view of (12. 23ft . Finally we pick p < r 
which gives the largest m. As t gets large, f](t,0) approaches —1 and p may be chosen so that 
m(t) = ming \ v {t, ?7(i, 0) + pe ie ) \ grows like t, namely, for p close to 1/2. In the limit t 00, fj(t,z) 
becomes holomorphic in the entire complex plane. 




2 



-2 



-4 



Figure 2: Profile of (12.141) for t = (solid line), 10 (long dashes), 10 5 (short dashes) and 10 20 (dots) 
and p in a neighborhood of p = —1 
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To describe the asymptotic behavior of r](t,z) as t — > oo, equation (12.151) can be written as 

"<'•«> = 2?-^ i {^ + ln( -' i) + 2i + ln ( 1 " i ? ie "' 



which gives 



i + 2t(r ? +l) + 0(t(r ? +l) 2 ,(r ? + l)) = z 



fj(t,z) = -i-^(±-z)+n(t.z) 



where, by Theorem 12.61 R is a regular function of z for |z| < m(£) which goes to faster than 1/t, 
concluding the proof of Theorem 12.21 Note that 



-x = \z\ 



,2 



i r 

lim lim —U(t,yNz) = lim w(t, x) = / \im p(t,x') dx' 

t—>oo N-^oo N t^oo J q t— >oo 

and £/o(z) = |z| 2 is an equilibrium solution of (11.121) . for any number of components N. 



□ 



Remark 2.7 Figure® shows the solution v(t,p) of \2.ll\) . for various t. Fort = 0, v(0,p) = t>o(p) 
is a monotone increasing (decreasing) function of p G (—4,0) (p G (—00,-4),) and z£s inverse 
v 1 (x) = u' (x) is defined for x G (—1/16, oo). Fort > 0, there is a unique negative value —l(t) (with 
1(0) = 4:), given byv p (t, —l(t)) = 0, such thatv(t,p) is monotone decreasing if — oo < p < —l(t) and 
monotone increasing if —l(t) < p < 0. The inverse function f _1 (t,x) = u x (t,x) has two branches 
but only the one with t> -1 (0,x) = u' Q (x) converges to —1 in any compact interval inside (—d(t), oo) 
with —d(t) = v(t, —l(t)) < —1/16 for allt > and d(t) — > oo as t — > oo. 



3 Geometry of the Scaling Flow 

Critical Trajectory. The scaling flow u(t, x), defined by equations (11.111) and (12. ip . is the cummu- 
lant generating function of the block spin variable at scale t. The flow is determined by its partial 
derivative u x (t,x) (see (12.101) ) and Theorem 12.21 exhibits a single trajectory, in the (viscosity) limit 
N -> oo, 

<ff(u' -> -1) = {u x (t, x), t > : u x (0, x) = u' Q (x), u x (oo, x) = -1} , 

that starts at t — from the initial function (12.131) and converges, as t goes to oo, to the stationary 
solution —1, implicitly defined by (12.151) and (12.91) . In this subsection we identify the class of 
functions where the flow is defined and give a geometric function theory description of this trajectory 
that establishes a one-to-one and onto relation between the orbit &(u' Q — > —1) and the time 
dependent convex domains Q(t) = u x (t,M), t > 0, formed by images under u x of the upper half- 
plane H = {z = x + iy G C : y > 0}. Analytical and numerical techniques are combined in order 
the conformal equivalence between Q(t) and H to be explicitly verified for all t. 

Analytic Continuation of initial value. Let us begin by extending Watanabe's proof of Propo- 
sition 12.11 to the upper half-plane HI. 
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Proof of Proposition HO . Let (j) u (£) = £,J U (£) / Jv-i(£) be defined for v > 1 and £ G C. The Bessel 
recursion relation 

J„_i(0 + ^+i(0 = 

generates a continued fraction of Gauss (see Chapter XVIII of |Wa] ) : 

2 (e/2) 2 2 (£/2f 

^ = ^ 1 = y 1 777^ ( 3 - X ) 

2u + 1 1 .J , m 



uniformly convergent over the domain 



1 lei 2 < i , (3.2) 



v(u + l) 

by Worpitzky's Theorem (see [Wa], p. 42). 

Let $ N {x) := U(0, VNz)/N with a; = - \z\ 2 . Equation (031) together with ([LTD and the Bessel 
recursion relation vJ v (£) — £J' V (£) = £J v+ i(£), gives 

1 f (N/2 - 1) Jjv/2-i (iyfflxN) - i^fi^N J' N/2 _ x {i^ffix~N) 



X§n{x) = 2N \ J N/2 ^{i^N) 

= ^<f>N/2(iy/Px~N) . (3.3) 



We take £ = iy/fixN and v = N/2 in (EH]) and write 

1 - r (f>N/2(i\^xN) 



a 2 



As iV goes to infinity, 

-fix 

ak = 



2k\ ( 2k + 2 
1 + — 1 + 



N J V N 

converges to — fix uniformly over the domain (13.21) for any integer k > and, consequently, x , d' N (x) 
converges over the same domain to a periodic continued fraction. We thus have 

xu>(x) = jim = ^ ^ = 1+ "f*^ (3-4) 

1 + -^- 
1+ '•• 

where the third equality is —1/2 times the solution of 

6 

Y 1 + 
14 



that is positive for positive x. This yields ( 12. 3 j) in view of the normalization $at(0) = 0. Note that 
the limit holds for any x in the domain 

\4{3x\ < 1 (3.5) 

of complex plane and this is sharp for the limit function uq since (— oo, —1/4/3] is a branching cut 
of u' . 

Lee— Yang zeroes. Both functions $' N (x) and u' Q (x) can be analytic continued to the upper half- 
plane and extended, by reflection, to the slit domain C\(— oo, —1/4/5]. As 0„(£) is a meromorphic 
(even) function of £, it can be written as 

^o = eY, a 2 1 _^ ( 3 - 6 ) 

n>l u i,f-l * 

where o; n)iy , n > 1, are zeroes of the Bessel function J v . So, the limit N — > oo of (13. 3p together with 
the asymptotic behavior of the Bessel's zeroes, 



a n , N/2 -i ~ (N - 1) - + (2n - 1) £ 



for n large, gives 



Wn(x) = lim fi'vix) = lim — rr ^ — = — / tt- — (3.7) 

° 2iV^^ 2tt J lf , -g(s) - (3x 1 ; 

for some positive function g satisfying g(s) ~ s 2 for large s. Note that {a n> jv/2-i, n > 1} are the 
Lee- Yang zeroes of the "a priori" initial measure (jl.7p and, by (11.131) and (11.151) . they become 
dense over an interval of real line. 

Pick class of functions. Let P denote the class of functions 

f(() = u(() + iv(0 , ( = x + iy, 

analytic in the upper half-plane EI with positive imaginary part: > if y > (see e.g. [Do] . 
Chap. II). The class of functions P forms a convex cone and is closed under composition: 

1. af x + bf 2 eP 

2. ho he P 

hold for any a, b > and h, h £ P- 

A linear function a + b(, a G R and 6 > 0, and the function — 1/C are clearly in P since both 
are one-to-one and onto maps of H into itself. It thus follows by ( 13. 61) . together with the properties 
1. and 2., that 4> u (() is in P and, in the topology of uniform convergence on compact subsets of HI, 
the sequence ($' N ) N>1 converges to i?^, in P ([Do], Sec. 4 in Chap. 77). Note the following equality 
^00 ( x ) = u o( x ) m the domain ( 13. 5ft and u' is the composition of four Pick functions: 1 + /3(, \/C, 
(1 + C)/7? and —1/C- This implies that u' £ P and equality between first and last expression in 
(13. 4|) holds with x replaced by £ e H, concluding the proof of Proposition 12.11 
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□ 



Integral representation. A function f(() = u (£) + iv (() is in the Pick class if and only if has a 
unique canonical integral representation [DoJ 



dfjt(X) (3.1 



where a = lim / (iy) jiy > 0, b = u(i) is real and \l is a positive Borel measure on R such that 



y^oo 



J (A 2 + l) 1 dfi(X) < oo. In addition, 



p ((a, b)) + = lim 1 f v(x + iy) <fe (3.9) 



2 J/iO 7T y a 

holds for any finite interval (a, b) and determines /i uniquely from /. 

The initial condition u' (x) of the flow u x (t,x) goes to as x goes to infinity (in any direction 
of the complex plane). Consequently, a of its canonical representation vanishes. In addition, b can 
be identified with the second integral. So, if 

/o(C) = jH J^d^W ( 3 - 10 ) 

is defined with dp(X) = p(X)d\ an absolutely continuous measure w.r.t. the Lebesgue measure dX: 

1 V4(-A)-l 

P(A) = 4^ (-A) (3 - U) 
whose support is — oo < A < —1/4, then 

u 'o( x ) = Phifix) 

/oo 
X^ 1 dfi(X) = —1/2 agrees with /o(0), by an explicit integration. Equation 
-oo 

(I3.10p together with (13.71) leads to the following relation between (13.111) and the empirical density 
\fg\s) of Lee- Yang zeroes: 27rp(A) = \(—g' o g^ 1 ) (A)] . 

Geometric function theory. By the Riemann mapping theorem (see e.g. [GK] ) if an open set 
Q is topologically equivalent to H (i. e. Q and EI are homeomorphic) then Q is also conformally 
equivalent to EI and there exist a biholomorphic (holomorphic one-to-one and onto) mapping / 
from Q to EL In some cases / can be made uniquely defined by Q. The conformal equivalence of 
open sets provides qualitative informations on the trajectory G{v!^ — ► — 1). As a function of the 
Pick class, u' maps EI into itself but we can be more specific about the image of H by m[,. From 
here on we fix (3 at the critical value /? c (4) = 4. We denote the upper semi-disc of radius r centered 
at xq by 

§ r ( Xo ) = {( = x + iy eM: (x - x ) 2 + y 2 < r 2 } 
and let 5^ t be the class in P indexed by t G R + satisfying 



16 



(i) if is an univalent function (one-to-one) 



(ii) f(Ct) — (t for some complex number ( t 

(iii) ^(1/2) = -1 

(iv) ^(c) = m 



Proposition 3.1 u' maps the upper half-plane H conformally into the interior of the upper semi- 
disc of radius 2 centered at —2: 

w' (H) =tt = S 2 (-2) 

and no other function in S^n maps HI into Qq. Hence, there is a one-to-one and onto relation 
between Q and the initial function u' Q of critical trajectory &{u'q — > —1) in the class S^q of functions 
with fixed point Co given by the complex root of2x 3 — x — 2. 



Proof. By equations (12.151) and (12.6j) . the inverse of u' , given by 



«b(p) = > ( 3 - 12 ) 



is the initial condition ( 12. 12ft of the linear evolution equation ( 12.111) . Hence 

Q = { v = p + iq e H : 9 (v (p + zp)) > 0} (3.13) 
and this is equivalent, by ( 13. 12}) . to the following inequalities 

q (p 2 - q 2 ) - 2pq (p + 2) > 0, q > 

which can be written as the upper semi-disc §^(— 2): {jp + 2) 2 + q 2 < 4, q > 0. 
Since 

«o(p) = 2p3~ 

does not vanish neither diverges for any p in Q but at edge points p = —4 and p = in the closure 
fio of fio, we conclude by ( 13. 12ft that tt (H) = Q is one-to-one and onto map. Note that u'q(x) 
vanishes at x — fo(poo) with p^ such that f (poo) = oo, in view of Mq o vo(p) = l/v' (p), i.e., at 
infinity in every direction of the complex plane. 

Now, suppose there exist another function f(x) in such that <£>(H) = $1q- Then, o u x 
is a map from HI onto itself, belongs to the class P and leave the points 1/2, Co and Co fixedQ As 
a consequence of (13.81) . the Pick functions that map H onto HI are linear fraction transformations. 
Since the identity mapping is the only linear fraction transformation leaving three points fixed, we 
infer that (p(x) and u x (x) are the same function. The complex root Co — —0.582687 + 0.720119i of 
the fixed point equation 

MO = -t^t = C 



4 The class of functions in P considered can be analytically continued across the real line by reflection (see 
condition (iv) of S*t)- If Co = ^o + Wo G HI is a fixed point of / € P then Co — — Wo is a fixed point of its 
extension. 
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is in Q , concluding the proof of Proposition 13. 11 



□ 

We now apply (13.131) to determine fl t = u x (t, H) for t > 0. As u x (t,() solves v(t,rj) = ( for 
rj, with v explicitly given by (I2.14p . the domain Q t can be easily plotted using ContourPlot or 
ImplicitPlot packages in Mathematica. Approximate expressions can be given for t around and 
oo. 

Domain boundary. Each set Q t of the family for t > is bounded by a simple convex closed 
curve which is piecewise analytic and defined by equation 

3 (v(t, r?)) = 0, 7] = p + iq e 1 (3.14) 

where v(t,r]) is analytically continued to the closure HI of the half-plane BL One has to be careful, 
however, in order to get the actual domain since Q (v(t,p + iq)) > may have more than one 
component. Figure [3] shows level curves of S (v(t,p + iq)). 



,A - 

l / 




Figure 3: S (v(t,p + iq)) = c for t = 0.2 with c taking negative (dashed lines), positive (solid 
curves) and neutral (thick solid line) values 

For t small, Qt is a slight deformation of §2( — 2), by continuity: 

/ 2(l + 2t)\ 2 9 4(1 + 2t) 2 
V 1 + 4t J (1 + At) 2 

Whereas, for t very large, Q t approaches a folium (half-leaf) of Decartes: 

3 f 2t l +P + Zg ) > 0, q > 2p (p 2 + q 2 ) + 3p 2 - q 2 < , q > . 
V (p + iqf J 
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We observe that the boundary of Q t is the union of two curves: a line segment I a := [—a, 0] 
extending from a point —a = —a(t) < up to the origin over the real line and a convex curve 
q = h(t,p) defined for p e I a with h(t,—a) = h(t,0) = 0. From the above, a(t) is a monotone 
decreasing function of t with a(0) = 4 and limf^ 00 a(t) = 3/2 whereas h(t,p) is a semi-circular 

curve at t = 0: h(0,p) = y 4 — (p + 2) 2 and approaches a limit (half-leaf) curve 



h*(p) = lim h(t,p) 



' 3p 2 + 2p 3 
1 - 2p 



Figure H] shows domain boundaries for various t. Note that Qt C fif if t' < t with strict inclusion 
along the convex arc. 




Figure 4: Domain boundaries Q t for t — n/4, n — 0, . . . , 9 

The function h(0,p) and the turning point — a(0) = —4 are related to the density p(X) of the 
canonical representation 113. lOft of u' (x) and its support E , 



P 



Q)=^V /Z ¥ :I 7=^(0, P ), -4< P <0 



by substituting A = 1/p in 1 13. lip . Note that S = E(0) = (— oo, — d(0)) in this case is such that 
d(0) = (4a(0)) _1 = 1/16. 

To determine the support E(i) = (— oo, — d(t)) of the measure p,(t,d\) = p(t,X)d\ of the 
canonical representation of u x (t,x) we look at the negative value —lit) at which v p {t, —l(t)) = 0. 
In the neighborhood of this point v is not univalent. Observe that l{t) and the turning point a(t) 
coincide. Writing 

v(t,r}) = y(t,P,q)+iw(t,P,q), V = P + *?, (3-15) 
by definition of a and Cauchy-Riemann equations, we have 

= w q (t, -a(t), 0) = y p (t, -a(t), 0) = v p (t, -a{t)) 
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which implies ot(t) = l(t) by uniqueness. From Remark 12.71 we have —d(t) = v(t, —l(t)) < —1/16 
for all t > and 

-d(t) ~ v(t,-3/2) 

9 27 V 2 / 27 

for t large enough, implies that the support £(i) of p,(t, A) converges to an empty set: £(t) = 
(— oo, —d(t)) — > as t — > oo. 

Riemann surfaces. Contour plots of v(t,rj), r] G C, for various t, show that u x (t,() is a mul- 
tivalued function of ( G C. Already at £ = 0, Uq(C) has two Riemann surfaces connected by a 
branch cut along the segment (—oo, —1/16] across which the imaginary part of u' Q (Q change sign 
(u' ((— oo, — 1/16]) is the semi-circular boundary of §2( — 2)). The determination of y/ ~ is chosen 
such that -1/ (1 + VI + 16C) 

is in P. For t > 0, has an even more elaborate Riemann 

surface with three sheets. The first is connected with the second sheet by a branch cut (— oo, — d(t)] 
while the latter is also connected to a third sheet by a branch cut (0, d\(t)\ with d\(0) = and 
d\(t) —>■ oo as t — > oo, which does not concern us as it doesn't relate to the limit function —1. The 
curves u x (t, (— oo, —d(t)]) and (0, di(t)]), which define together with the real line boundaries 
of two domains, intercept the real line perpendicularly at negative and positive values, respectively. 
Figure [5] shows these curves for various t. The region bounded by u x (t, (0, di(t)]) inside the half- 
plane EI is denoted by A t . Note that, opposed to Q t , A t are open domains satisfying inclusions 
A t c A t ' if t < t'. 




-4 -3 -2-1 12 



Figure 5: Domain boundaries VL t and A t for t = n/4, n = 0, ...,9 

Flow in the Pick class. It is very difficult to show directly from the flow equation that u x (t, C) 
remains in the Pick class of functions for all t > by general principles. However, for initial 
condition in P that belongs to the class J?^ there is a simple property of the flow equation that 
explains why the Pick class P is preserved. Writing vq(t]) = yo(p,q) + iw (p,q) as a function of 
rj = p + iq G fi U Qq U I a (o), with Q* the reflection of Q about the real axis, if the imaginary 
part Wq(p, q) is an odd function of q then the flow equation (12. lip preserves this property. Writing 
v(t,T]) as (I3.15p . we have 

w(t,p,q) = -w(t,p, -q) 
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holds for all t > and i] £ Q t - By continuity, it follows that Q t G M and that v(t,rj) remains a 
one-to-one and onto map from Q t to HI and these imply that u x (t, Q belongs to the class in P. 




Figure 6: Scaling evolution of fixed points and domain images of u x (t, Q 

To establish uniqueness of the relationship between 0(u' Q — > —1) and the image domains 
{Q t , t > 0}, we proceed as in t = (see proof of Proposition 13. ip . Supposing that <p(t, £) G S" t 
is a different function satisfying (p(t, H) = Q t , for each t fixed (p~ l o u x (t, () maps H into itself 
and leaves the point 1/2, ( t and £t* fixed where Q = ( t for t < t* ~ 5.155075 and for t > t* the 
last two fixed points become real numbers (see Figure [6]). Extending the functions in across 
the real line by reflection, o u x (t, () is a linear fraction map with three fixed points which 
contradicts the hypothesis that u x (t, () and (p(t, () are different. This holds for all t such that 
Ct G fit U O* U I a {t)- If this condition is not satisfied, we apply Schwarzian reflection |Daj about 
the curve h(p) in order to extend u x (t, () to the complex plane in such way that u x (t, H) = Q t and 
u x (t, — H) = H\(fi t U At) and this insures that ( t and Q, which are now real values, remain fixed 
points of u x (t, () when Q < —a(t). The value t co that Q co = —a(t co ) is called crossover scale from 
strong to weak (coupling) regime, term introduced in [HHWj. 

The canonical representation of u x (t, () is not suitable for describing the trajectory <ff{u' Q — > — 1). 
From the characteristic equations (12.161) of (12. lip one find that the point (p, V) = (—1,1/2) is a 
critical point for the two-dimensional dynamical system: 

(p, \>) = (-2p (1 +p) , -1 + (6 + Ap)V) := (F 1 (p, V), F 2 (p, V)) (3.16) 

with Fi(-l,l/2) = F 2 (-l,l/2) = 0. As (-1,1/2) is an invariant point we have v{t, -1) = 1/2 
and, accordingly, u x (t, 1/2) = —1. Instead of fixing b in the canonical representation (I3.8P the 
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value of / at £ = i, we write 

«. (*,0 = -1 + £ - ^) «W«, A) ■ (3.17) 

Note that, by hypothesis, [(A — £) (A — 1/2)] _1 d/i(t, A) < oo and as the support £(t) = 
(— oo, — of fi(t,X) converges to the integral in (I3.17P converges to uniformly in each 
compact set OeH. 

The following summarizes our findings. 



Proposition 3.2 u x (t, (), t > 0, map the upper half-plane H conformally into a decreasing family 
of open convex sets Sit satisfying 

Sl t = u x (i,H) c «o(e) = 

and no other function in S?t maps H into Sl t . There is a one-to-one and onto relation between this 
family and the trajectory &{u' Q — > —1) at the critical inverse temperature (3 = /3 C (4) = 4. The geo- 
metric description together with the integral representation ofu x (t,x) gives the distribution d/i(t, A) 
of the Lee-Yang zeroes at the scale t. Sloo is a nonempty set and a nontrivial limit distribution is 
attained but its support £(£) is pushed away from the origin to infinity. 



4 Normal Fluctuations 

We turn our attention to normal fluctuations. The block variable (11.81) is now normalized with 
7 = d and the system is above the critical temperature. The "a priori" measure (x) , L dK = n, 
that governs the law of (jl.8p . satisfies a recursive equation 

a[ N \x) = -L e ^(^i)MV2 V £0 * ... * afJ x {L^x) , k > 1 

S v ' 

L d — times 

which, in view of 7 = d, has an explicitly dependence in the exponential pre-factor (see (II. 3p ). 

Initial value problem. Following the procedure described in Section the initial value problem 
(I1.12p and (I1.13p . for the logarithmic of its characteristic function 4>^\z) in the L J, 1 limit, thus 
reads 

U t = - l -e- 2t (AU - \U Z \ 2 ) +dU- 1 -z-U z + ^e- 2t AU(t,0) . (4.1) 

Note that, L~ 2k = exp (— 2A;lnL) — ► exp (— 2t), as k — > 00 together with L I 1 with fclnL = t fixed, 
and such function appears in front of the Laplacean in (11.61) . 

As N — > 00, the radially symmetric solution of (14. ip scaled properly satisfies the modified initial 
value problem (see (11.151) for the definition of u(t, x)): 

u t = e~ 2t u x - 2xe~ 2t u 2 x - ^xu x + du - e~ 2t u x {t, 0) (4.2) 

with u(0, x) = Uq(x) given by (12.31) . 
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We continue through equations (I2.5I) - (I2.10I) . A similar Legendre transform applied to (14. 2 p leads 
to the initial value problem 

v t - 2p 2 e- 2t v p = -e~ 2t +(d + 4e" 2 V) v (4.3) 

with v (0,p) = v (p) as given by (I2.12p . Note the cancellation of terms proportional to pv p because 
7 = d in this case. 

Main result. The following result holds for any d > 2 by it has been stated for d — 4, for 
simplicity. 

Theorem 4.1 Equations jjj4.3\ ) and Ii2.12\) with d = 4 are solved by 



_ e 4t 



v (t,p) = -^-(l-^ + l-ln(l- P + pe~ 2t ) + e~ 2t (in (l - p + pe' 21 ) - l) - V 4 ^ . (4.4) 



For every f3 < /3 C (4) = 4 and t > 0, there is a unique solution p = p(t, x) of 

v(t,p) = x , (4.5) 

holomorphic in a neighborhood of the origin, that converges exponentially fast, as t —>■ oo, to the 
solution of 

4 p 

in every compact set ofC This implies, together with the corresponding equations Ii2.10\) and \2. 
convergence to a Gaussian equilibrium solution of the equation ( [^.1| ) without terms proportional to 
e- 2t : 

lim lim \-U (t, VNz) = (4.6) 
t^ooN-+ooN V / 2/i(/?) v ' 

uniformly in compact subsets of — \z\f G C. 

Proof. As in the proof of Theorem 12.21 equation (14. 3 p will be solved along the characteristics 
p{t;po). We refer to this proof for details. Writing V{t) = v(t,p(t)), we have 

p = -2e~ 2t p 2 

V = -e~ 2t +(d + Ae- 2t p) V (4.7) 
with initial conditions (12.171) . Integrating the first of these equations gives 

1 + Po - Poe 2t 

The homogeneous equation V = (d + 4e~ 2t p) V can be integrated analogously as before 

V(t) = V e dt (l+p -p e- 2t ) 2 
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Using the variation of constants formula, the solution to the second equation of (14.71) is given by 

V(t) = e dt (l +Po - Po e~ 2t ) 2 (V - Jo) (4.9) 

with 

1 f 1 C d/2 d( 
J o — : 



'exp(-2t) (1 +Po - PoC) 

At this point, notice that Jo for d — 4 is exactly as in the proof of Theorem 12.21 Equations (12.171) . 
(12. 121) together with the integration of Jo gives 

V + J = (JZ2Q]) . 

The difference between the two cases is the exponential pre-factor e dt of (14. 9p and p = po(t,p) 
which is now obtained by solving (14.81) for p : 

Po (t,p) = - V —— Tt . (4.10) 

1 — p + pe zt 

As we shall see, these two differences are responsible for the converge of trajectories to different 
stationary solutions. 

Equation (14.41) follows by plugging (12.2011 into ( 14.91) with p given by (14.10j) . 

We now solve equation (12.151) for p at (3 ^ (3 C = 4 which, by (12.14[) . can be written as 

(xp 2 - |) e~ 4i - e~ 2t = -1 + £ - (V 2t - ±\ In (l - p + pe^) = 9l (t,p) . (4.11) 



Analogously to Lemma I2.5[ we have 



Lemma 4.2 For any p < (1 — e 2< ) , g\ is a monotone increasing function of p with gi(t, 0) = 
— e~ 2t +/3/4 and diverges logarithmically to — oo as p — > — oo. 



Proof of lemma. g\ is a monotone increasing function of p since 

P -2t fi _ p -2n i 

fa),( t .") = 1 _ P ( 1 _ <; J) + ^/(i»(i -«-*)) 

with / given by (12.221) is a positive function for .p < (1 — e~ 2t ) 1 . Other statements follows as in 
the proof of Lemma 12.51 

□ 

The quadratic polynomial Qi(x,p) in the left hand side of (14. lip tends to a linear function 
Q!(0,p) = -e" 2 *(l +pe~ 2 V2) as x -> with Qi(0,0) = -e" 2 * and Qi(0, -2e 2 *) = 0. From Lemma 
14.21 the graph of g\ always intercepts the graph of Qi(0,p) for all /3 > and, as in the proof 
of Theorem 12.21 this implies the existence of a unique solution p(t,x) of ( 14. 5ft for every t > 0, 
holomorphic in a neighborhood £7(t) of the origin that becomes the entire complex plane U(t) — > C 
as t — > oo. Details of the proof will be omitted since are similar to the corresponding statements 
in Theorem 12.21 
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Asymptotic expansion. The asymptotic behavior of p (t, x) as t — > oo is given as follows. By 
equation (14. lip . p(t,x) converges exponentially fast 

to a constant value p which solves 

l-£ = zl In (1-0 = ^(0) . (4.12) 
4 p 

Since /ii is a monotone increasing function of p < 1 with /ii(0) = 1 and lim^oo hi(p) = 0, there is 
a unique solution for all < (3 < 4. Comparing (I4.12p with (jl.lOj) . together with (11.111) . f)2.ip and 
( l2~TUj) . equation gSJ holds with 

1 

P ~ MP) 

concluding the proof of Theorem 14.11 

□ 



5 Conclusions and Final Remarks 

In the present work, a continuous version of the hierarchical spherical model at dimension d = 4 
has been investigated. The two main results are Theorems 12.21 and 14.11 on the limit distribution of 
the block spin variable X 1 normalized with exponent 7 = d + 2 at the criticality and 7 = d above 
the critical temperature. To prove these results, certain evolution equations corresponding to the 
renormalization group transformation (11.61) in the limit L | 1 are solved explicitly at iV = 00. 
Starting far away from the stationary Gaussian fixed point the trajectories of these dynamical 
system pass through two different regimes with distinguishable crossover behavior. The large-iV 
limit of the transformation (II .6p with L d fixed equal to 2, at the criticality, has been investigated in 
both weak and strong (coupling) regimes by Watanabe [W] . We mention that our analysis using the 
L [ 1 limit equation is considerably simpler and, consequently, has more details than Proposition 
2.2 in |Wj. 

Theorem 13.21 gives an interpretation for the above mentioned trajectories using the geometric 
function theory. The methods used enable us to describe the dynamics of the Lee- Yang zeroes 
along those trajectories. As iV — > 00, the Lee- Yang zeroes becomes dense over a semi-line and 
their measure, which depends on the scale parameter t, is shown to reach a limit for t large but the 
support of the limit measure is pushed away to infinity as the trajectories approach the Gaussian 
fixed point. The method also allow us to give the precise crossover scale t co > t* from strong to 
weak regime defined as the value of t such that Q = —a(t) where (* is a fixed point of the function 
(I2.14p that solves equation (12.111) and — a(t) is a point of the boundary of image domain Q t . 

There are, however, two major drawbacks in the L j 1 limit equation of the hierarchical O(N) 
Heisenberg model with iV finite. Firstly, reflection positivity cannot be used to prove uniform 
convergence of the O(N) trajectories to 0(oo) trajectories. 

The other problem is related with the Lee- Yang property. A Borel measure p in M> N possesses 

Lee- Yang property if its characteristic function <p(z) — I dp (x) exp (iz ■ x) belongs to the Laguerre 
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class J?f of entire function of ( = — |z| 2 G C which can be represented by 

/ (C) = exp (AC) (l + ^) (5.1) 

oo 

with A > and ai, a 2 , • • • real numbers satisfying a^ 2 < oo. Hence (see [Nj IHHWt IW] ) 

fe=i 

oo 

MO = "C (In /)' (NC) = £ (-l) j C J (5.2) 

is holomorphic function of £ in a neighborhood of the origin and Newman's inequalities 

< u 2j < (is 4 y /2 (5.3) 

holds for all j > 2. The scaling (15.21) is chosen so that v 2 j — O (1) in N for j > 1 if p is the uniform 
measure on the sphere of radius yN. Inequalities ( 15. 3p can be shown to hold in the limit AT — > oo 
but in this case / cannot be represented by (15. ip as the zeroes (ctj) >1 become dense over the real 
line. 

Now, let 

fk = Tfk^i, k — 1, 2, . . . 

where T : £ — > <§ is the operator defined by recursion relation (II. 6p with f(Q = ip (\z\) = (f>(z), 
( = — |2;| 2 , be a sequence in the space of entire functions § starting from /o(/3C) with / in the 
Laguerre's class «5f. It has been proven in Theorem 1.1 of |KWj that, for every k G N and 
< < (L^- d - 1)/A, 

fk g =$f n srf\ 

where srf a denotes the Frechet space of functions / G $ such that 

1 



keN o k 



d k f 
} (0) 



k 



is finite for all b > a and A is the type of fa. This together with equation (15. 3p can be used to 
establish the existence of a critical inverse temperature (3 C such that the sequence {/ n } ngN converges 
to exp (() uniformly in compact subsets of C. The Pick class of functions is the natural candidate 
for replacing Laguerre's class in the local potential approximation of (11.61) but we don't have a 
substitute for the convex space s$\. The present work is an attempt in this direction for N = oo. 
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